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Abstract: This paper deals with the attitude estimation of a rigid body equipped with angular 
velocity sensors and reference vector sensors. A quaternion-based nonlinear observer is 
proposed in order to fuse all information sources and to obtain an accurate estimation 
of the attitude. It is shown that the observer error dynamics can be separated into two 
passive subsystems connected in "feedback". Then, this property is used to show that 
the error dynamics is input- to- state stable when the measurement disturbance is seen as 
an input and the error as the state. These results allow one to affirm that the observer 
is "robustly stable". The proposed observer is evaluated in real-time with the design and 
implementation of an Attitude and Heading Reference System (AHRS) based on low-cost 
MEMS (Micro-Electro-Mechanical Systems) Inertial Measure Unit (IMU) and magnetic 
sensors and a 16-bit microcontroller. The resulting estimates are compared with a high 
precision motion system to demonstrate its performance. 
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1. Introduction 

1.1. Motivations and Background 

The attitude control problem of a rigid body has attracted a strong interest during the past few decades, 
and it has been extensively studied. This interest comes from the fact that many aerospace systems, 
such as spacecrafts, satellites, tactical missiles and many others, enter within the framework of a rigid 
body with the requirement of precise attitude information [1]. In the last decade, the application of 
micro-electro-mechanical systems (MEMS) has gained a strong interest. Therefore, the attitude 
estimation problem has been tackled within new areas, such as terrestrial and aerial robotics [2,3], 
virtual reality [4], biomechanics [5] and wearable robots [6]. In these cases, the attitude information 
is obtained from inertial and magnetic sensors, namely, three rate gyros, three accelerometers and three 
magnetometers, orthogonally mounted, such that the sensor frame axes coincide with the principal axes 
of the rigid body. Since the attitude is not directly measured, it must be estimated via the measurements 
of the mentioned sensors. In general, these sensors can be classified into two main categories: Angular 
velocity sensors that measure the angular velocity of the body with respect to some inertial frame and 
reference vector sensors that give the coordinates of a fixed vector in the mobile frame. 

Rate gyros are angular velocity sensors, and they provide continuous attitude information with good 
short-term stability when their measurements are integrated by means of the rigid body kinematic 
Equation. However, since rate gyro measurements are affected by drifts and noise, the attitude estimation 
based on these sensors diverges slowly from the real attitude. On the other hand, when reference vector 
sensors measurements, known as vector observations, are available, the attitude can be determinated by 
algebraic or optimization techniques [7]. Therefore, in order to decrease the effect of noise induced by 
the vector observations, as well as the attitude divergence provoked by the rate gyro drifts, it is important 
to fuse both information sources. 

To this effect, various estimators have been proposed, using the quaternion attitude 
parametrization [8]. The Extended Kalman Filter (EKF) [9] and new alternatives to the standard 
EKF have been applied extensively (see [10] and the references therein). However, the divergence 
problem [11], the non-Gaussian noise induced [12] and the computational cost render difficult the 
embedded implementation of the Extended Kalman Filter (EKF). In recent years, significant research 
efforts have been addressed toward the synthesis of nonlinear attitude observers in order to tackle 
the previously mentioned issues. The first work on this topic was presented in [13]. Subsequently, 
in [14-16], quaternion-based nonlinear observers and gyro bias observers are proposed in the framework 
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of satellite sensors calibration and marine vehicle navigation, respectively. Furthermore, the formulation 
of nonlinear attitude observers based on the the orthogonal group, SO (3), have been proposed in the few 
last years [17-20]. Recently, an excellent overview of rigid body attitude estimation and comparative 
study was given in [21,22]. Nevertheless, it is well known that a nonlinear observer is generally 
vulnerable to measurement disturbance, in the sense that a small arbitrary disturbance can lead to the 
divergence of the error state [23]. The above mentioned nonlinear attitude observers (with the exception 
of [20]) do not consider either the present noise in the gyro measurements nor a time- varying gyro 
bias term. Furthermore, the computational complexity in these approaches is a drawback, since the 
algorithms work with 3x3 matrices and they have to preserve orthogonality properties. The singular 
value decomposition (SVD) is a well-known approach and an extremely powerful and useful tool in 
linear control theory and signal processing. However, few works have exploited this approach in the 
framework of attitude estimation. In [24], the authors use the SVD method to estimate the attitude matrix 
minimizing Wahba's cost function. In the work reported in [25], the authors use the SVD approach for 
finding the relative orientation of a robotic manipulator. Nevertheless, the problem of attitude estimation 
from magnetic and inertial sensors using an SVD approach has not been addressed in the literature. 
Furthermore, unlike sensors and computer systems used in the marine and aerospace community, the 
signal output of the low-cost sensors is subject to high levels of noise and a time- varying bias term. 
These problems must be addressed, since, in a practical framework, the design of efficient and embedded 
attitude estimator using low-cost sensors is an important issue. 

7.2. Contributions 

In the present work, a quaternion-based attitude observer/estimator of a rigid body is presented. In the 
proposed approach, the attitude estimation problem is solved in two parts. Firstly, a quaternion attitude is 
estimated by means of vector observations. In this first step, the attitude estimation is performed using an 
SVD (singular value decomposition) approach. Then, the quaternion obtained in this step is considered 
an attitude measurement. Contrary to conventional techniques, the SVD maintains the quaternion's unit 
constraint naturally. Furthermore, the numerical robustness and numerical stability are guaranteed [26]. 
The second part of the proposed method consists of the design of a nonlinear observer in order to produce 
an estimate of the time-varying gyro bias and the attitude quaternion. This observer is driven by an 
attitude error obtained by means of the quaternion propagated by the observer, and this one obtained 
from the SVD technique. Asymptotic convergence of the estimation error is proven. Moreover, it is 
shown that the error dynamics can be decomposed in two passive subsystems connected in "feedback". 
This result is exploited to prove that the observer is input-to-state stable (ISS) [27,28] when the rate gyro 
noise is seen as the input and the error as the state. In this sense, using the small gain theorem, one claims 
that the observer is "robustly stable". To evaluate the proposed attitude observer behavior in real-time, 
a complete Attitude and Heading Reference System (AHRS) based on low-cost inertial and magnetic 
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sensors and a 16-bit microcontroller is designed and implemented. A comparison with a high precision 
motion system is carried out, in order to demonstrate the observer performance. 

The ISS paradigm in an attitude observer, the problem of estimating the attitude from vector 
observations using an SVD approach, as well as a real-time implementation have never been addressed 
in the literature. These facts show the originality of the present work. 

The document is organized as follows. In Section 2, a mathematical background of the attitude 
parametrization and sensors modeling is given. The main problem is formulated in Section 3. The 
formulation of the nonlinear attitude observer and stability analysis is presented in Sections 4-6. The 
AHRS implementation and experimental results are given in Section 7. Finally, conclusions and further 
research are mentioned in Section 8. 



2. Mathematical Background 

2.7. Unit Quaternions and Attitude Kinematics 

Consider two orthogonal, right-handed coordinate frames: the body coordinate frame, 
E 6 = [e$, e^, e|], located at the center of mass of the rigid body, and the inertial coordinate frame, 
= [e{, ef, ef], located at some point in the space. The rotation of the body frame, E 6 , with respect 
to the fixed frame, E f , is represented by the attitude matrix R G 50(3) = {R G M 3x3 : R T R = 7, 
det7?= 1}. 

The cross product between two vectors, x,£ G M 3 , is represented by a matrix multiplication 
[£ x ]x = <f x Z where: 

( o -6 6 

[e x ]= 6 o -a 

V -6 6 o J 

The n-dimensional unit sphere embedded in W n+1 is denoted as § n = {x G : x T x = 1}. 

Members of 50(3) are often parametrized in terms of a rotation, /3 G M, about a fixed axis, e G § 2 , by 
the map, Z7 : R x 8 2 -> 50(3), defined as: 

«(& e) := 7 + sin(/?)[e x ] + (1 - cos(/?))[e x ] 2 (1) 

Hence, a unit quaternion, q G S 3 , is defined as: 



2 

e sin | y \ q 



»==! ) = [* 1^' (2) 



<7 = [?i ?2 Qs] T € K 3 and g 0 £ R are known as the vector and scalar parts of the quaternion, respectively. 
q represents an element of 50(3) through the map, 1Z : S 3 — » 50(3), defined as: 

n:=I + 2q 0 [q x ]+2[e x } 2 (3) 
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Note that R = 1Z(q) = IZ(-q) for each q e S 3 , i.e., quaternions g and — q represent the same physical 
attitude. The inverse unit quaternion is given by q~ x := [q Q — ^ r ] T , and the quaternion product is 
defined by: 

( Ql ° ~^ \ ( Q2 ° \ fA\ 

gi®92 := I T _ X1 ^ (4) 

Denoting by uj = [uji uj 2 oos] T the angular velocity vector of the the body coordinate frame, E 6 , 
relative to the inertial coordinate frame, E^, expressed in E 6 , the kinematics equation is given by: 

= ~ f \o = \mo <5) 
a) H h% + m ) 2 

The attitude error is used to quantify the mismatch between two attitudes. If q defines the true attitude 
quaternion and q is the estimated quaternion, then the error quaternion that represents the attitude error 
between the true orientation and the estimated one is given by [8]: 

q e := q~ l <g> q = (q eo q e T ) T (6) 

In the case that the true quaternion and the estimated one coincide, the quaternion error becomes: 

q e = (±1 0 T ) T (7) 

2.2. Sensor Modeling 

The sensors in the rigid body are classified into the two following categories. 

2.2.1. Angular Velocity Sensors 

The angular velocity uo = (uj\ cj 2 u 3 ) T is measured by three rate gyros orthogonally mounted. The 
output signal of a rate gyro is influenced by noise and by a slowly varying function, v, called bias [29], 
that is: 

uj g = uj + V + ff G (8) 
V = -T- l v + fj y (9) 

where ffc and fj u G M 3 are bounded noises and T = r/ 3 is a diagonal matrix of time constants. 

2.2.2. Reference Vector Sensors 

Let us consider the representation of a vector, x u with respect to E^ and E 6 and denoted by fi and 
bi, respectively. The vectors, f i9 are also called the "reference vectors" and, in general, are known quite 
accurately. The body vectors, bi, are known as "vector observations" and are obtained from on-board 
sensors (for instance, magnetometers, star trackers, accelerometers). Let b iq and r iq be the associate 
quaternion to the vectors, 6- and ?v 

K = (o K T ) T n q = (o n T ) T (io) 
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these quaternions are related by the quaternion rotation, q, such that: 

K =q~ 1 ®r lq ®q (11) 
From Equation (4) and using quaternion algebra, the following measure model can be obtained [25]: 

0 ~(b i -r j ) T 



q = Hq = 0 

with H G M 4x4 . Note that Equation (12) is a well- structured linear system of equations. 



(12) 



3. Problem Statement 

The objective is to design an attitude observer that estimates the rigid body attitude doing a trade-off 
between a good short-term precision given by rate gyro integration and a reliable long-term accuracy 
provided by vector observations. 

If at least two vector observations are non-collinear, the quaternion can be estimated using these 
simultaneous vector measurements. For this purpose, the measurement Equation (12) is exploited. The 
observation matrix in a function of m pair vectors (6 i? fty with i G {1, 2, m} is written as follows: 



/ 



H 



\H m ) 



o -(h-n) 
h-n -[(6i + fi) 



0 



V 



( bfn ) 

■[(b m + f m ) x ] 



p4mx4 



(13) 



/ 



In the case of perfect measurements, the measure Equation Hq = 0 is always satisfied. However, 
in practice, the measurements are polluted with perturbations originated from noise, vibrations, etc. 
Therefore, the problem consists of estimating the attitude quaternion from a perturbed observation 
matrix, H, that is: 



(14) 



q ps = argmm<j ^ \\Rq[ 2 
subject to \\q ps \\2 = 1 

On the other hand, if the initial value, q(to), is known and the angular velocity measurements 
are perfect, i.e., Co = ujq, it is possible to obtain q(t) by integrating the kinematics Equation (5). 
Nevertheless, q(t 0 ) is, in general, unknown, furthermore, the angular velocities are polluted with bias 
and noise (see Equation (8)), which results in slow divergence of attitude over time. Therefore, the 
observation idea can be used to on-line correct the integration, q(t), of some erroneous q(t 0 ) and attitude 
errors caused by the bias, and the noise presents in the rate gyro measurements. This correction can be 
carried out according to the measurable error q e = q~ l ® q ps , where q ps is considered an attitude measure 
obtained by the vector observations. 
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The observation problem can then be formulated as follows: Given the kinematic equation of the rigid 
body and the model of time- varying bias evolution, find an estimate q(t) and v{t) from the knowledge 

of uo{t) and q ps (t) for all t >t 0 G M>o. 

4. Nonlinear Observer Formulation 

4.1. Attitude Estimation from Vector Observations 

— * — # 

Assumption 1. There are at least two linearly independent vectors, bi and bj, with i ^ j, that is, 
rank(H) > 3. 

Lemma 1. Let H G ]R 4mx4 be the observation matrix formed by the m pair of vectors (^,f^) with 
i G {l,2,...,m}. Then, there is an orthonormal matrix U = {U\ U2 ... U± m ) G ]R 4mx4m , an 
orthonormal matrix V — [V\ V2 V3 V4) G M 4x4 a/irf a^z orthonormal matrix S G R 4mx4 , wzY/i elements 
ah along the diagonal and zeros everywhere else, where k — 1 : 4, ^wc/i ^tftf: 

^ = [/^ T (15) 

Then, the attitude quaternion is given by the last column of the matrix, V, that is: 

q ps = v 4 (16) 

furthermore, the normality condition \\q ps \\2 1 is satisfied. 

Proof. The proof follows the standard least squares solution of homogeneous equations via the 
SDV approach [30]. □ 

4.2. Nonlinear Attitude Observer with Bias Estimation 

The proposed attitude nonlinear observer that includes the bias and the error update is given by: 

[ V = -T~ X V - K 2 q e 

where T G M 3x3 has been defined in Equation (9) and K\, K 2 are positive constant parameters, q is the 
prediction of the attitude at time t. It is obtained via the integration of the kinematics Equation (17), 
using the measured angular velocity, ujq, the bias estimate, v and q e , which is the vector part of the 
quaternion error, q e , that measures the discrepancy between q and the measured attitude, q ps , obtained 
by means of Lemma 1, that is: 

q e = q~ X ® q ps = (q eo fef ( 18 ) 

Actually, the observer Equation (17) allows one to fuse the data that arise from the vector observations 
and the rate gyro measurements. The observer structure is shown in Figure 1. 
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Figure 1. Block diagram of the nonlinear observer (® quaternion multiplication, 
© matrix multiplication.) 
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Note that the attitude quaternion, q ps , can be expressed by: 

q ps = q® Aq 



(19) 



where q denotes the true attitude and Aq represents a perturbation quaternion, which is caused by the 
sensor measurement noise encountered in practice. 

For the first time, a Lyapunov stability analysis is carried out. As is always done, the error dynamic 
is obtained in a noise-free context; then, the following hypotheses are made: 

• Vg — Vv — 0- The measurement rate gyro noise is not included in the error dynamics. 

• Aq = (1 0 0 0) T , i.e., q = q ps . This is a good assumption for the case of static attitudes and for 
quasistatic movements [17]. 

Combining Equations (5)-(9) and (17), the error dynamics is expressed as: 



0 



7 



S e : { 



-7 [2^ X ] + [7 X ] 




(20) 



T~ x v + K 2 q e 



where 7 = v + K\q e and v = V — V. 

The system Equation (20) admits a set of equilibrium points V := {p e +, Pe-}, where: 

• p e+ := (1 0 T 0 T ) T 6 § 3 x R 3 



p e . := (-10 T 0 T ) T G§ : 



3 x R 3 



In both scenarios, one shows that the following stability property is ensured with respect to 
different sets. 
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Definition 1. (Asymptotic Stability in the Large) [31]. Let X C M n be given. The trivial solution x = 0 
ofx = /(t, x) w called asymptotically stable in the large with respect to X if it is stable in the sense of 
Lyapunov and every other solution x(t, to? x o) ~^ Oast —> oc for any initial states, x 0 G X, and for any 
initial times, £ 0 £ K>o- 

One considers that the state space of Equation (20) corresponds to either 

X+ = {q e G § 3 | q eo G [0,1]} x M 3 or X_ = {q e G § 3 | q eo G [-1,0]} x M 3 , depending on 
which target equilibrium is chosen. 

Proposition 1. Assume that all trajectories with initial conditions, (g eo (t 0 ) Q^(to) v T (t 0 )) T G X+, 
satisfy sgn(q eo (t)) = sgn(q eo (to)) > 0 for all t > t 0 (mutatis mutandis for X_). Then, the equilibrium 
point, p e+ (respectively, the equilibrium point, p e _) of the system Equation (20) is asymptotically stable 
in the large with respect to X+ (respectively, with respect to X_). 

Proof. Consider the candidate Lyapunov function, V : S 3 x IR 3 —> R, which is proper and 
positive definite: 

V = K 2 ((l - q eo ) 2 + f e q e ) + h T v = 2K 2 (1 - qj + \v T v (21) 
The derivative of Equation (21), together with Equation (20), is given by: 

V = ~2K 2 q eo + v T b = -K 2 fq e + z> T (-T- x z> + K 2 q e ) ^ 
= -K 2 {f + K r q T e )q e - v T T~ x v + K 2 v T q e = -K 2 K x qJ q e - v T T~ x v < 0 

Thus, q e —> 0 and v —> 0 and, due to the normality condition, q eo —> 1. Hence, solutions of the system 
Equation (20) whose initial conditions are in X+ converge asymptotically to p e + := (l 0 T 0 T ) T . That 
concludes the proof. □ 

Remark 1. In a practical context, one aims at improving the convergence speed of the observer. Since, 
physically, the set of equilibrium points, V := {p e+ ,p e -}, corresponds to the same attitude error, it 
is necessary to establish the desired equilibrium point to be achieved, depending on the given initial 
condition. This improvement can be ensured by choosing the equilibrium point corresponding to the 
sign of q eo (to). The point, p e +, is chosen if q eo (to) > 0 and the point, p e _, otherwise. Then, a slight 
modification on the observer Equation (17) is done, obtaining 



0 Adis :{ : 2 S( "V~ -™ -> (23 ) 



Q = 2 S W) ~ ^ + sign(q eo )K 1 q ( }j 
v = -T~ x v - sign(q eo )K 2 q e 



This fact can be shown adapting the proof of Proposition 1 using the following Lyapunov 
function instead: 

v = . K 2 ((l ~ q eo ? + £ Qe) + \v T v, if q eo (t 0 ) > 0 
K 2 {(1 + q eo ) 2 + gq e ) + \v T v, if q eo (t 0 ) < 0 
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5. Passivity Interpretation of the Attitude Nonlinear Observer 



Lemma 2. (Kalman-Yakubovich-Popov (KYP) lemma) Let Z(s) = C{sl — A)~ X B be a p x p transfer 
function matrix, where A is Hurwitz, (A B) is controllable and (A, C) is observable. Then, Z(s) is 
strictly positive real (SPR) if and only if there exist a positive definite symmetric matrix P = P T > 0 
and Q = Q T > 0, such that: 



A T P + PA = -Q 
B T P = C 



(25) 
(26) 



Since K 2 is a scalar, the dynamics error Equation (20) for the bias observer can be written as: 

b = -T~ l i> + K 2 hq e 



(27) 



Although the error system Equation (27) has no output, it is possible to define the following 
virtual output: 

y e = K 2 hv (28) 
Using Equation (28), the error dynamics Equation (20) can be structured as the one shown in Figure 2. 

Figure 2. Block diagram of the dynamics of attitude and bias error. 



o 



Hi 



—T~ 



I3K2 



Proposition 2. The transfer function matrix of the error system Equations (27) and (28) is SPR. That is, 
the mapping, q e — >> y e , satisfies the KYP lemma. 

Proof. Since the matrix, — T -1 , is symmetric and a Hurwitz, one has for all Q symmetric positive definite 
matrices a P symmetric positive definite matrix, such that: 

{-T~ l ) T P + P(-T~ l ) = -Q (29) 
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If the Q matrix is specified as Q = 2T 1 , one gets: 

P = / 3 (30) 

then, 

£ T P = {K 2 h) T h = K 2 h = C (31) 
which satisfies the equality Equation (26). □ 

Proposition 3. The dynamics of the quaternion error, represented by the system, %\, is state 
strictly passive. 

Proof. Consider the candidate Lyapunov function, V q , which is positive, definite and proper: 

v = [K 2 {{\ - q eo ) 2 + <f q e ) = 2K 2 (1 - qj, si q eo > 0 
(K 2 ((l + qj 2 + £q e ) = 2K 2 (1 + qj, si q eo <0 

Analyzing for q eo > 0, the derivative of Equation (32) along the trajectories of q e (Equation (20)), is 
given by: 



V q = -2K 2 q eo = -K 2 fq e 



= -K^v 1 + K x qi )q e = -K 2 K x q e L q e - K 2 v L q ( 

or, similarly: 



(33) 



K 2 v T q e = V q + K 2 K x qJq e (34) 



Since q eo e [0, 1], one obtains: 



q e T q e = l-ql>l- q eo > (1 - q eQ f (35) 



ie 0 — ^e 0 — V- 1 - ^e 0 

Then: 

^2^§ > + ^K^l -q eo f (36) 

Power flow Rate of change of State dissipation 

energy stored rate 

□ 

Theorem 1. The attitude nonlinear observer Equation (17) is passive. 

Proof. Propositions 2 and 3 establish that the systems, H 2 and Hi, are SPR and state strictly passive, 
respectively. Since the two systems are connected in feedback (see Figure 2), the attitude nonlinear 
observer Equation (17) is passive. □ 
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6. Input to State Stability (ISS) of the Attitude Nonlinear Observer 

Note that the convergence and the passivity interpretation of the observer were accomplished within 
a context free from noise and perturbation. However, the rate gyro measurements are corrupted by 
measurements noise and bias components that are time- varying. Owing to the nonlinearity of the 
observer, this type of disturbance can cause instability or even the finite time escape of the estimate. 

In this section, the robustness of the proposed observer is discussed, which follows from passivity 
propriety. It will be shown that the error dynamics is input- to- state stable when the disturbances are 
viewed as the input and the error quaternion as the state. 

In the case that noise is included, the error dynamics Equation (20) becomes: 



1 



0 



Qe 



7 +Vg 




(37) 



2 1 -7 + ?7g [2£ x ] + [(l + ffoY 
i) = -T~ x v + K 2 q e + rj v 
This equation can be described as illustrated in Figure 3, where the inputs, d\ and d 2 , are such that: 

di = v + fj G , d 2 = K 2 q e + r\v 

It is assumed that the noise components are bounded. That is: 

(38) 

However, neither the bounds nor the noise distribution are necessarily known. 



sup ||77g(s)||2 < ci sup ||^(s)|| 2 < c 2 

se[t 0 ,t] se[t 0 ,t] 



Figure 3. Block diagram of the dynamics of attitude and bias error with 
measurement disturbances. 



VG 



T(7 + jfc,^) 



H, 



Qe 0 



:0 



hK 2 



Proposition 4. The system, % 2 , is input-to-state stable (ISS) with respect to the input, d 2 . Then, the 
output ofH 2 (y e = K 2 I 3 i>) is always bounded. 
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Proof. Consider the second equation of the system Equation (37). Since the matrix, — T 1 ,isa Hurwitz, 
one obtains: 

= o-r- 1 (*-*o)iirw 



\Hto)\\2+ [ e- T -^-^d 2 (s)ds 
= e-r-^Mto)^ + fe- T -^- s \K 2 q e (s) + vAs))ds (39) 

< e -T-( t - t o)||^ o) || 2 +R2 sup {||gr( a )|| 2 }+ SUP {U(S)\\2} 

se[t 0 ,t] se[t 0 ,t] 



then: 



\\v(t)\\ 2 < ma,x{( 2 (\\i>(t 0 )\\ 2 ,t-t 0 ), sup <r 2 i(||<f e (s)|| 2 ) , sup <? 2 (||^(s)|| 2 )} (40) 

se[to,t] se[t 0 ,t] 

where: ((\\i>(t 0 )\\ 2: t - *o) = e -r-i(*-* 0 ) | |i>(* 0 ) 1 1 2 is a KC function and <r 2 i(||<f e (s)|| 2 ) = ^ 2 ||<7 e (s)|| 2 and 
s" 2 (||?tL(s)|| 2 ) = ||?7i/(s)|| 2 are /C functions. Then, the bias estimation error, v{t), is input- to- state stable, 
and the output y e (y e = K 2 l^v) is always bounded. □ 

Proposition 5. 77?e system, Hi, is input-to-state stable with respect to d 2 . 

The key to showing the ISS of the system, Hi, is the fact that a system is input- to- state stable if it 
admits an ISS-Lyapunov function. In this sense, a Lyapunov function will be proposed; then, it will be 
shown that this function is a ISS-Lyapunov function. 

Proof. Consider the candidate Lyapunov function, V, which is proper and positive definite: 

(K 2 ((l-q eo f + gq e ) = 2K 2 (l-q eo ), si q eo > 0 

Vg = < (41) 

(K 2 ((l + q eo ) 2 + £q e ) = 2K 2 (l + q eo ), si q eo < 0 

Analyzing for q eo > 0, the derivative of Equation (32) along the trajectories of q e (Equation (37)) is 
given by: 

V q = -2K 2 q eo 

= -K 2 (f + ffJ)q e 

= -K 2 {y T + K x f e + ffj)q e (42) 
= -K 2 K x qJq e - K 2 {v T + ffj)q e 
= -K 2 K x qJq e - K 2 d^q e 

since q eo G [0, 1], one obtains: 

Qe T Qe = MeWl = 1 ~ > 1 ~ <7e 0 > (1 ~ Qeof (43) 

then: 

V q <-K 2 Ki{l-q eo f-d^q e 

< -^l(l-(?eo) 2 +l|^l||2 ||<fe|| 2 

<-K 2 Ki{l-q eo f + \\di\\ 2 

< -a 3 ((l - q eo f) + a4(|Mi|| 2 ) 



(44) 
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since V q is proper 3 a l5 a 2 /Coo functions, such that: 

«i(lk(to)|| 2 )<K(^o))<« 2 (||<?(to)|| 2 ) 

Furthermore, a 3 ((l — g eo ) 2 ) = KiK 2 ((l — q eo ) 2 and ce^ ( 1 1 c?i 1 1 2 ) = ^IMilb are /C^ functions. Then, 
the Lyapunov function, V q , is a ISS-Lyapunov function (see [27]), and this fact implies that the system, 
Hi, is ISS. □ 

Since the system, Hi, is ISS, the following inequality is accomplished: 



\qe(t)\\ 2 < max<{ Ci(l|9e(*o)ll2,*-*o)> SU P Sd(\\di(s)\\ 2 ) 

se[to,t] 

< max<{ Ci(l|9e(*o)||2,*-*o) 3 sup <a 2 (||i/(s)|| 2 ), sup Si{\\ff G (s)\\ 2 ) 

s€[to,*] s€[to,*] 



(45) 



where d is a /C£ function and <*i 2 (||z>|| 2 ) = ^Gll^lh and ^1 ( 1 1 rfo 1 1 2 ) = ^ 2 ||?7g|| 2 are JC functions. That 
is, every trajectory of the attitude error dynamics asymptotically approaches a ball around the origin, 
whose radius is a function of supremum norm of the bias estimation error and of the measure noise of 
the rate gyro. 

Using the interesting properties of ISS systems, the stability of the entire error dynamics is studied. 
Then, one has the following result. 

Theorem 2. Let K 2 be a positive constant, such that K 2 < 1. Then, the attitude nonlinear observer is 
robust towards to measurement disturbance of the rate gyro. 

Proof. The claim follows from the small gain theorem for the interconnection of nonlinear systems that 
are input- to-state stable [32]. The error dynamics of a nonlinear observer (Equation (37)) is composed 
of the feedback interconnection of the systems (Hi-H 2 ). From Propositions 4 and 5, it follows that: 

|9e(*)||2 < max^ Ci(l|9e(*o)||2,*-*o), sup & 2 (||i?(s)|| 2 ), sup ^(\\vg(s)\\ 2 )\ (46) 

se[t 0 ,t] se[to,t] 

||z>(t)|| 2 <max<jC2(||z>(t 0 )||2,t-to), sup s 21 (\\q e (s)\\ 2 ) , sup s 2 (\\ff„(s)\\ 2 ) \ (47) 
[ se[t 0 ,t] s£[t 0 ,t] J 

The composition of the functions, s 12 and is given by: 

a2 0fti(||ge(s)|| 2 ) = K 2 (K 2 \\q e \\ 2 ) (48) 

Since K 2 < 1, it follows that: 

K 2 \\Qeh< Meh V ||<f e || 2 >0 (49) 

Then, the system (Hi-H 2 ) is ISS, when the measurements disturbances, fjc and ff u , are seen as the 
input and the error state, \qj £ T ] T , as the state. This means that the attitude observer is robust with 
respect to measurement disturbance. □ 
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7. Experimental Results 

The estimation methodology proposed in this work is implemented and evaluated in real time, in order 
to assess its effectiveness. For this purpose, an embedded system was designed and developed. Special 
attention was paid to the low power consumption requirements and weight, leading to the selection 
of the digital signal controller (DSC), dsPIC33FJ128MC802, which was used with a clock speed of 
4 MHz. It contains extensive digital signal processor (DSP) functionality with a high performance, 
16-bit microcontroller (MCU) architecture, but without a floating point unit. The sensor suite is based 
on a sensor board equipped with a tri-axis accelerometer (ADXL135), a dual axis gyro (LPR530AL), 
a single axis gyro (LY530ALH) and a tri-axis magnetometer (Micromag 3). All sensors outputs are 
analog, except for the Micromag 3, which is digital and uses the serial peripheral interface (SPI) bus 
system as the underlying physical communication layer. The aim was to test the DSC implementation 
with an update rate of the estimation of at least 55 Hz. 

Furthermore, the system is equipped with a Bluetooth module (BlueSMiRF Silver), which provides 
wireless capabilities. The total system supply voltage is 3.3 V. The dimension and weight are 
60 x 40 x 15 mm and 60 g, respectively. The system is depicted in Figures 4 and 5. 

Figure 4. Block diagram for the Attitude and Heading Reference System (AHRS) prototype. 
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In order to experimentally evaluate the performance of the proposed AHRS, the experiments were 
carried out with a two degrees-of-freedom high accuracy system. Actually, this system is the Twin Rotor 
System [33] (see Figure 6). This system has two rotational joints, which allow for movement about two 
orthogonally axes with a resolution of 0.35° considered a true attitude for performance evaluation. In 
order to test the complete attitude estimation, two trials were carried out. In the first one, the axes, and 
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e 3 6 , of the AHRS coordinate frame (body coordinate frame E 6 = [e$, e|, e|]) coincide with the axes, 
and eg , of the Twin Rotor System. Thus, 6 (pitch) and ip (yaw) angles were evaluated. In order to evaluate 
the angle, cj) (roll), the AHRS was rotated tt/2 about the axis, e\ 6 , such that e\ 6 coincides with and e 2 6 
with e 3 * . Note that the attitude is estimated with respect to the inertial coordinate frame = [e{, ef , ef ] . 
E^ is chosen to be the NED (north, east, down) coordinate frame. In this case, the reference vectors 
are the gravitational and magnetic vectors. The vector observations, i.e., the gravitational and magnetic 
vectors expressed in the body frame, E 6 , are obtained from the mentioned tri-axis accelerometer and 
tri-axis magnetometer. The angular velocity is obtained from three rate gyros. 

Figure 5. The AHRS prototype. 





Figure 6. Experimental setup. 




The attitude and gyro bias observer is implemented at 55 Hz using a fourth-order Runge-Kutta 
method, with tuning parameters, Ki = 3.2 and K 2 = 0.9, such that the restriction of Theorem 2 is 
accomplished. No particular emphasis was given on the tuning process, since the resulting performance 
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with these simple parameters is very acceptable. The attitude is maintained as a quaternion, which is 
normalized at each time step to compensate for numerical errors. After that, the quaternion is converted 
to Euler Angles, since they are more intuitive, and it allows for comparing with the angles obtained from 
the Twin Rotor System. For the entire experimentation time span, the vector iff = v Q + V r was added to 
the vector of the rate gyro measurements, in order to test the ability to estimate and compensate for the 
rate gyro bias. In this case, one has V Q = 0.1(1 — 11) T and V r = 0.0001fc(l 1 1) T , where k represents 
the iteration index. 

Figure 7. Estimation of 6 (pitch) and (yaw) angles. 
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Figure 8. Estimation of 0 (roll) and i/j (yaw) angles. 
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For the first trial, the initial condition for the quaternion estimate is q(t 0 ) = (0.85 0.06 0.11 — 0.49) T , 
i.e., (0 = 0°, 9 = 15°, $ = —60°). For the second trial, the initial condition for the quaternion estimate 
is q(t 0 ) = (0.76 0 0 - 0.64) T , i.e., (0 = 0°, 6 = 0°, %j) = -80°). In both cases, the rate gyro bias 
initial estimate was set to zero. The observer angle estimates are compared with the ones obtained from 
the encoders of the Twin Rotor System. In Figure 7, the evolution of 6 (pitch) and x/j (yaw) angles is 
depicted. Figure 8 shows the evolution of the second trial, i.e., (f) (roll) and (yaw) angles. Finally, the 
convergence of the rate gyros bias estimates is shown in Figure 9. These results indicate that the observer 
faces up to the slowly time- varying nature of the rate gyro biases. 

Figure 9. Evolution of the rate gyro bias estimate. 
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It is worth mentioning that the algorithm itself uses only 81% of the 55 Hz cycle time. 
This is depicted in Figure 10. Finally, we invite the readers to watch the following video: 
https://www.dropbox.eom/s/r3pq57rwlf0pxuv/AttitudeObserver.mov. 



Figure 10. Algorithm CPU usage. 
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8. Conclusions 

This paper presented an original quaternion-based strategy for the attitude estimation of a rigid body 
equipped with angular velocity sensors and reference vector sensors. To accomplish the objective, a 
robust nonlinear observer for attitude estimation was developed. This observer ensures the asymptotic 
estimation convergence and local exponential estimation convergence. Furthermore, by means of an ISS 
analysis, it was shown that the observer is robust towards noise and the disturbance of the rate gyro 
sensors. Thus, an embedded AHRS was designed and developed. Experimental results were presented, 
where the AHRS outputs are compared with a two degrees-of-freedom high accuracy motion system. 
This illustrates the achievable performance of the proposed estimation schema. 
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